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Abstract 



We explore the four dimensional Einstein-Maxwell landscape as a toy model in 
which we can formulate a sphere compactification stabilized by an electromagnetic 
field. Replacing the compactification sphere by J spheres, we obtain a simple sector of 
the (2 J + 2)-dimensional Einstein-Maxwell landscape. In this toy model, we analyze 
some properties which are very difficult to uncover in the string theory landscape, 
^-j. ■ including: complete moduli stabilization, stability conditions, and state counting. 

A detailed comparison between the main features of this landscape and the Bousso- 
Polchinski landscape is given. We finally speculate on the impact of these phenomena 
■ in the string theory landscape. 



1 Introduction 



As a candidate of a theory-of-everything, string theory has led to many striking results. 
Among them we find the string theory landscape [1] [2] , a very complicated structure of 
vacuum states of the theory which raises its own questions and problems. All models of 
this landscape are rich and complex, and the existence of this landscape is almost beyond 
doubt |3J. When deriving cosmological models in a landscape we are led to the notion 
of multiverse [3], a quantum ensemble of different classical cosmological models. The 
multiverse shows very appealing features, such as inflation, which is needed to remedy the 
difficulties of older cosmological models [5j El [7J El |9] , and it is generally believed that the 
cosmological constant [TOj E] and the coincidence problems can be solved with realistic 
models of the multiverse [T2l[T3] . Nevertheless, the multiverse has its own problems, which 
will be described briefly. 
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The huge amount of possible universes present in the landscape [HJ UHl UHl E] should 
be complemented with a probability distribution which can explain why we are living in 
this particular universe. The standard rule of asigning probabilities as proportional to e~ s 
for some euclidean action S breaks down in this context because the classical action of 
a cosmological model is divergent due to the infinite spacetime volume. To extract some 
useful information, the action should be regularized, but there are many different ways of 
regulating an infinite volume because no cut-off procedure is invariant under a coordinate 
change. Different regularization procedures can lead to different probability distributions 
predicting different universes. This is known as the measure problem [H [18] , which consists 
in giving an unambiguous definition of the relative probabilities in a given landscape model. 
So far, the measures derived from first principles need the AdS/CFT correspondence [19] 
and similar ideas as a key ingredient [20j [21], which has an stringy origin, and thus the 
whole picture is slowly evolving towards a unified formalism of the multiverse derived from 
string theory. 

Another problem is the vast complexity of the landscape. Extracting a four- dimensional 
cosmology from string theory requires choosing a compactification of the remaining dimen- 
sions. Each possible choice of compactification describes a sector of the entire landscape, 
and the number and properties of the states in different sectors can be very different. This 
diversity can be explored by considering different sectors separately. Thus, simplified mod- 
els of some sectors of the landscape have been built, and they constitute a very important 
tool to understand the full implications of the string theory landscape. The reader can 
consult [22j [23] for a review of several models of the string theory landscape. For our 
purposes only the simplest models are needed. Among these models we should mention: 

• The Kachru-Kallosh-Linde-Trivedi (KKLT) mechanism [2J1 |25l [26] is a supersymmet- 
ric string model with many stabilized AdS vacua which are lifted to dS by quantum 
effects. This lanscape model gives inflationary states represented by brane-antibrane 
pairs evolving in Klebanov-Strassler throats in the compactification manifold. The 
model has no classical dS vacua, and thus its more controversial aspect is the quantum 
nature of the lifting from AdS to dS, which is uncontrollable by its very nature. 

• The Bousso-Polchinski (BP) model [1] is a simplified vacuum sector of M theory 
compactified from eleven to four dimensions with quantized four-form fluxes and 
M5 branes wrapped around three-cycles in the seven- dimensional compactification 
manifold. The volumes of the three-cycles are the moduli of the model, and gives 
a fundamental charge which determines the cosmological constant of the 3+1 
cosmological part of the eleven- dimensional metric by the formula 



Here, each is the number of flux units stored in the i three-cycle, and A is 
a bare, negative cosmological constant, a parameter of the model. Equation ([I]) 
provides a very elegant mechanism to deal with the cosmological constant problem. 





i=l 
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By choosing a large number J of three-cycles, we can choose the integers rii to 
approximately cancel the A contribution, thus obtaining a huge landscape of AdS 
and dS vacua containing states with very low values of the effective cosmological 
constant A. Calabi-Yau compactifications typically have a large amount of three- 
cycles, so that realistic values of the cosmological constant can be obtained without 
any fine-tuning of the model parameters. 

On one hand, unlike the KKLT mechanism, the BP model does not rely on quantum effects 
for producing dS states. On the other hand, BP moduli are frozen from the very beginning, 
and thus no stabilization mechanism analogous to KKLT is included in the BP model. 

A very appealing feature of the BP landscape is that the closed expression (JTJ) leads 
to tractable counting problems. Thus, relative probabilities of different states can be 
computed based on combinatorics, but this leads to the conclusion that the number of 
states with realistic A values are a very sparse minority, thus opening the door to anthropic 
arguments J27J [23 12] . This is a consequence of the measure problem in the BP multiverse, 
that is, of the lacking of a preferred and completely specified prescription for computing 
probabilities. 

If we combine the computation of probabilities by naive state counting with the lack 
of stability analysis, we find another curious property of the BP landscape. Under a 
large number of three-cycles J, flux space is a very high-dimensional space. Such spaces 
have very counterintuitive geometric features: for example, the opposite vertices of a unit 
hypercube are separated by a distance yJ, which can be interpreted (in Planck units) as 
a huge energy scale. Another consequence of this elongation of diagonal distances is that 
the vast majority of lattice nodes inside a sphere in flux space are located on hyperplanes, 
that is, they always have at least one vanishing component. In fact, we have seen that as J 
grows, the dimension of the most populated hyperplanes is distributed in a narrow gaussian 
window around a typical value a* J with a* < 1 [29]. But as the KKLT mechanism shows, 
flux quantum numbers cannot vanish in stable states. Of course, the lack of a stabilization 
mechanism in the BP model prevents us from directly excluding those states from the 
model; if they were excluded, the states with very low A values could even disappear. This 
"a*-problem" of the BP model is only present for large J, but a large J is certainly needed 
to solve the cosmological constant problem. 

Thus, looking for a scenario where we can closely examine this a*-problem, we are led 
to consider a model with the following properties: 

• It should be exactly solvable. 

• It should have many moduli. 

• It should have a stabilization mechanism. 

The simplest landscape with the first and third requirements is the Einstein-Maxwell land- 
scape, considered in a plethora of papers as a toy model landscape [3U1 ED E21 S3]- This 
landscape has compactifications of the form (A)dS 2 x S 4 or (A)dS 4 x S 2 having a single 



3 



modulus, namely, the radius of the compactification sphere. We have added many moduli 
by considering compact ifications of the form 



(A)dS 2 xnS 2 , 



(2) 



and thus the moduli are the radii of the J spheres. We call this sector multi-sphere 
Einstein-Maxwell compactification. We will see that dS states with low quantum numbers 
are always unstable, and thus they should be excluded from the landscape. This legitimates 
the a*-problem as an objection against naive counting arguments in the BP landscape. 

We have chosen the cosmological part in the multi-sphere model to be 1+1 instead 
of 3+1 for simplicity. On one hand, these cosmologies are unrealistic, and they have the 
peculiarities of two-dimensional dilatonic gravities. On the other hand, this multiverse 
is safe against anthropic reasoning, and thus a good candidate for studiyng the measure 
problem in future papers. 

The paper is organized as follows. In section |2} we give a detailed description of the 
single- modulus compactification (A)dS 2 x S 2 in the Einstein-Maxwell theory, together with 
the stabilization mechanism and the state counting. In section [3] we consider the multi- 
sphere model and its stabilization mechanism in detail, and give a very detailed account 
of the (A)dS 2 x S 2 x S 2 sector of the six- dimensional Einstein-Maxwell theory. Section H] 
is devoted to counting the states in the model. In section we summarize the differences 
found between the multi-sphere Einstein-Maxwell sector and the Bousso-Polchinski land- 
scape, and in section [6] we speculate on some consequences that those phenomena can have 
on the string theory landscape. The last section [7] summarizes our conclusions. 

2 One-flux compactification in the four-dimensional 
Einstein- Maxwell theory 

2.1 The one-flux four-dimensional Einstein Maxwell landscape 

By a compactification in four- dimensional gravity we understand a solution of the Einstein 
field equations of the form 



This is a particular expression of a Kantowski-Sachs cosmology [Ml 135] . The metric splits 
in a (t, x)-spacetime part, which we will identify with a two-dimensional cosmological 
solution, and a (u, w)-surface /C, the compact part. 

The Einstein field equations are, in units with G = c = 1, 



ds 2 = e 2H*,t) (_ dt 2 + ^2) + e ty M ^2 + ^ 



(3) 




(4) 
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The stress-energy tensor has two contributions, the first comes from the electromagnetic 
field and the second from a vacuum energy density: 



T {M) = _L( F p P _ I F 2 \ 

^ V 4tt\ w v 4 M V 



(5) 



T (A) = _ 



A 

8tt 



9nv 



The symbol F 2 = F^F^ is the electromagnetic Lagrangian density, and A is the cosmo- 
logical constant of the four-dimensional theory. 

The stress-energy tensor of the Maxwell field is traceless, so we have 

A 



T = T (A) = 



2tt 



(6) 



We will assume a magnetic monopole configuration for the electromagnetic field: 

F - ^(■-"'"'■''M„ A dr. (7) 



^e 2 ^duAdv, 



where the boldface is used to denote differential forms. Here the constant Q is the magnetic 
charge of the monopole, and V is the volume of the compactification manifold, so that 



and we have 



V = vo\lC= [ e 2 ^ u ' v) duAdv 
F = Q. 



(8) 
(9) 



K 



In matrix notation, the Maxwell field is 
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B x = — e 
V 



(10) 



This configuration solves Maxwell equations in curved spacetime 

V„F^ = -j= d u (y^F^ = 
and the non- vanishing components of the corresponding stress-energy tensor are 



(11) 



(M) 



— Bl e 2 ^ = — ( Q ) e 2(f) 



1 

8n Lx 



1 

8n\V 



1 



r(M) - D2 20-4^ 

8tt x 



U9l 



8n\V 



D 2<f> 



(12) 



rp[M) _ rp{M) 



(9l 



8n\V 



3 2V 



The contribution of the cosmological constant is 



Finally, the Ricci tensor has the following nonzero components: 

Rtt — —Rxx = —<t>tt + 4>xx , 

Ruu Ryv ^uu Ipvv 1 



(13) 



(14) 



where the subscripts in <p xx etc. represent partial derivatives. Einstein equations coincide 
in the tt and xx components, and in the uu and vv components also: 



~<Ptt + 4>xx = 



-A + 



2-1 



(15) 



Thus, and ip are uncoupled and satisfy Liouville equations of — h and ++ signatures 
respectively. 

The Liouville equation states that the Gaussian curvature of the corresponding surface 
is constant. We will call these two constants A and K: 



(0«-0„)e- 2 ^ = A-^ 2 = A, 



A + 



(16) 



Thus, a solution of the — h Liouville equation cf>(t, x) represents a two-dimensional space- 
time of constant curvature A, which is de Sitter (dS 2 ) if A > 0, Minkowski (M 2 ) if A = 
and anti-de Sitter (AdS 2 ) if A < 0. Therefore, A can be interpreted as the cosmological 
constant of the dimensionally reduced cosmological model. On the other hand, ip repre- 
sents a compact surface of constant curvature K. We can choose K positive, and then the 
surface will be a sphere of radius 1 / \[K. 

The positivity of the constant K is equivalent to zero genus by the Gauss-Bonnet 
formula, which also relates volume with curvature: 



— / Ke 2 ^dudv = 2 
2vr 



K 



KV n TA 4vr 
= 2 V = — , 

2tt K ' 



(17) 



which leads to an algebraic equation for K: 



Q 



K = A + ( ^ ] = A 



QK 
An 



(18) 
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The previous equation has two solutions: 

2 



*T ± = 2A' Q 



Q 



i± Wi 



Q 2 



Qr 



Qmax r-r- • (1^) 

v A 



The magnetic charge should not exceed Q max . For greater charges curvatures become 
complex and therefore the solution makes no sense. 

The two-dimensional cosmological constant has also two branches 

A± = 2A - K T . (20) 

Assuming the usual Dirac quantization condition on the magnetic charge in terms of the 
elemental charge e of the particles coupled to the electromagnetic field, 

Qe = 2im, withneZ, (21) 

we have a maximum value of the integer n, 



(22) 



and all integers n satisfying < \n\ < n max , plus n — 0, constitute the one-flux four 
dimensional Einstein-Maxwell landscape^. 

The case n = deserves further comment. In this case we have Q = 0, and we have 
only one branch, K — A = A, which cannot be supported because the electromagnetic field 
vanishes. Thus we should expect this solution to be unstable, as we will see below. 

Thus far, we have considered A > 0. In the case A < 0, it can be seen that the solutions 
of equation ( TT8|) yield K_ < and thus only the K + branch remains as a solution, with a 
cosmological constant A_ = 2A — K + always negative. In addition we have no restriction 
in the quantum number n, and therefore this A < infinite landscape is less interesting 
than its A > counterpart, which we will be considering henceforth. 

Equation ( 120]) can be interpreted as a distribution of the dS4 curvature A between 
the (A)dS 2 and S 2 parts. The solutions obtained show that a positive curvature A can 
be distributed between positive K and positive or negative A (two possible ways), but 
a negative curvature A should be distributed between a positive K and a negative A (a 
unique way), yielding a physically less interesting landscape. 

Figure [1] shows curvatures and cosmological constants for n max = 50 for both branches. 
Note that this number controls the number of the states in the landscape, whereas A 
controls the magnitude of the moduli K and A. 



2.2 Modulus stabilization 

The landscape considered thus far has only one modulus, namely the Gaussian curvature 
of the compact sphere or equivalently its radius. This is a volume modulus, which control 

1 This is the "pedestrian" landscape mentioned in Footnote 2 of Ref. [36] . 
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Figure 1: Curvature X (left) and two-dimensional cosmological constant A (right) of the 
two branches of the one-flux four-dimensional Einstein-Maxwell landscape are shown with 
"■max = 50, corresponding to A = 3.915 x 10~ 4 (a random value, expressed in units such 
that e = 1). The state n = is present only in the A+ branch. The branching point is 
absent because Q max is (generically) not an integer. The n = 1 point on the K + , A_ branch 
produces finite values K + (l) ~ — A_(l) ~ 10 4 A. The horizontal line signals the onset of 
stability, all states above it in the A+ branch being unstable. Thus, this landscape contains 
50 AdS2 states and 3 dS2 states. 



the volume (surface area) of the compact part of the spacetime. The next question we 
should ask is the stability of the solutions encountered. If a value of K in one of the two 
branches for fixed quantum number n is not stable, a small perturbation will drive the 
system out of the solution. If the subsequent evolution makes the sphere radius to grow 
unbounded we speak of decompactification. An unstable state in the landscape should not 
be included when counting states in the resulting two-dimensional multiverse. 

The one-flux compactified solution found is dS 2 x S 2 or AdS 2 x S 2 with a sphere radius 
which is homogeneous throughout the two-dimensional spacetime. The perturbation will 
alter this situation, assuming that the compactification radius can be different on different 
t, x points, yielding an ansatz 

The perturbation £(£, x) appears also on the non-compact part of the metric, and in this way 
the local volume element remains invariant. This reflects that the perturbation exchanges 
locally volume between the compact and non-compact parts of the metric. Functions <j) 
and if) are solutions of the Liouville equations encountered before. 

We will investigate the dynamics of the perturbation field £ (t, x) from a two-dimensional 
perspective. The four- dimensional action of the Einstein-Maxwell system is 

= _L / L W dt dx du dv , L (4) = y^^(R^ -2A- F 2 ) , (24) 
lo7r J 

where we use superscripts to distinguish between four- and two-dimensional quantities. 
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The curvature scalar and the electromagnetic scalar are 



J2« 



25! e" 4 ^- 4 « . 



(25) 



Using a/— g( 4 ) = e 2<?!>+2 ^ and integrating by parts the second derivatives of £, the four- 
dimensional Lagrangian is 



L (4) 



2£ 



+ (£-£)e* + ^ + Ae 2 * +2 * + B=e : 



20-2i/>-4£ 



(26) 



Now, we will substitute the sphere ansatz — (^i uu + tp vv )e~ 2 ^ = K and the magnetic 
monopole configuration .R,. = ^e 2 ^ as backgrounds for the dynamics of the perturbation. 
Integrating with respect to u, v we obtain the two-dimensional action 



S^=V- I L^dtdx 

4:71 



= 1 e 2 ^( 2 ) - I ^(g^Up - U(0 



(27) 



where the summation in a, (3 e {t,x} is implied, and 

(y 2) ) Q/3 da d/3 = e 2 ^" da d/3 = e 2 *(-dt 2 + dx 2 ) , 



u(0 



A-Ke~^ + 



QK 
4tt 



(28) 



that is, a model of 1 + 1 spacetime with gravity non-minimally coupled with a scalar (called 
the dilaton or the radion) which experiences self-interaction through a potential. Note that 
the dependence with respect to the dilaton is not relegated to the potential, and so a direct 
stability analysis using £/(£) is not possible. 

The next step is to show that the Euler-Lagrange equations of the previous two- 
dimensional model are equivalent to the Einstein equations of its four- dimensional coun- 
terpart, that is, the truncation of the model is consistent. Firstly, we recast the Lagrangian 
( 127]) displaying explicitly all fields and removing the second derivatives of by integrating 
by parts, which amounts to the substitution 



for a G {t,x}. The resulting Lagrangian is 



L (2) 



- e ^T] 



(29) 



(30) 



The equations of motion are 

V° P (^ - Cap - = -e 2 V(0 (31) 

with respect to £, and 

t P + e 2€ = -2e 2 ^[/(0 (32) 

with respect to (p. The absence of perturbation £ = should be a solution of the equations, 
so that equation (j3~T]) with £ = reduces to 

- e- 2 V^ = + &*) e ~ 2 * = ^'(O) = A > ( 33 ) 
which is the two-dimensional cosmological solution, and equation (1321) gives 

C/(0) = A-iT+(^) 2 = 0, (34) 



which is the same equation previously found for K, see (JT8J). We will now assume that 
£ is a small perturbation, thereby neglecting the backreaction of the perturbation on the 
geometry of the cosmological solution. This allows us to fix <fi as another background field 
by means of equation (1331) . This eliminates as a dynamical variable in the Lagrangian 
( |30l) and therefore equation ( l32l) will not be used. In other words, we assume that equation 
( |32l) is satisfied to zeroth order, which is the content of (134p . and we are left with ( 13T1) as 
the evolution equation for the perturbation £ in the background 0. 

Thus, we substitute eq. ( 133]) into eq. ( )3T|) . resulting in a dynamical equation for £ which 

is 

- e-V'feu + = A - e~ 2 ^'(0 = -C4(0 . (35) 

The linear stability analysis of equation (1351) requires its linearization (the effect of neglect- 
ing the non-linear derivative term does not spoil linear stability, as discussed in appendix 

ED 

e-^fet-U) =-U2f(«)e. (36) 
which is a 1 + 1 Klein-Gordon equation. We also consider a small spacetime region such 
that can be treated approximately as a constant. We thus obtain the linear stability 
condition of the solution £ = as being a minimum of the effective potential: 

t^(0)=4(21f-3A)>0, (37) 

which is K > | A. All points in the K + branch satisfy this condition, and therefore all 
AdS2 states are stable, but this is not so in the K_ branch. The condition K_ > |A is 
met by all charges satisfying 

2^/2 

Q > Qmin = g Qmax • (38) 



Upon quantization, the previous condition is 

Qr 



2tt 



(39) 
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Thus, all states in the dS2 branch characterized by a quantum number n < n m - m are 
unstable. This includes also the state n = 0, as mentioned above. We conclude that the 
flux number labeling true vacuum states should obbey a double inequality, obtained by 
combining (1221) and ( )39l) : 

^min < n < n max . (40) 
Therefore, the number of stable states in this one-flux landscape is 



AdS 2 



+ n Vl 



- n n 

dS 2 



+ 1 



e 




\2y/2 e i 






3 x/A 



+ 1 



2 - 



2^ 



(41) 



In the example shown in Figured] we have Q max = 317.54 (corresponding to A = 3.915 x 
10~ 4 ) which gives M\ = 53 (50 AdS 2 states and only 3 dS 2 states). 

The physical description of these states is simple: The self-gravitating electromagnetic 
field of the monopole stabilizes a geometry (A)dS 2 x S whose natural behaviour is decom- 
pactificatiorH, that is, the curvature of the sphere part, which tends to vanish, is sustained 
by the magnetic field. The distribution of curvature contributions between the compact 
and non-compact parts of the geometry is whatever allowed by magnetic charge quan- 
tization in the AdS 2 case, while in the dS 2 case only large charges can sustain positive 
curvature of de Sitter states. 



3 Adding many fluxes 

Obtaining a non-trivial landscape with many fluxes is not easy. The easiest technique is 
to extend the electromagnetic field F to a SO (J)-invariant J-component field, in which 
all charges are different [TJ [36] . This approach is not convenient to address the problem of 
stabilization, because the charges do not come from a transparent compactification scheme, 
and therefore nothing is known about the stabilizing potential. 

In all known compactifications, the charges come from moduli describing the shape of 
the inner manifold. These moduli are free geometric parameters, but they are promoted to 
dynamical scalar fields in the perturbation analysis. The charges are considered coupling 
constants, and therefore their dynamics should be governed by a potential with at least one 
minimum. The stabilization problem consists of finding this potential. Different models 
provide a wide variety of potentials; if the potential does not possess any minimum, then 
the dynamics of the moduli will lead them to grow unbounded; this phenomenon is known 
as decompactification. 

In an ideal model, we should expect that all moduli come from a compactification 
manifold which is derived from the equations of motion. Nevertheless, these equations are 
very difficult to solve in its full generality, and therefore the inner manifold is chosen at 

2 Note that equation (|36|) . when the state is unstable, as happens in absence of electromagnetic field, 
predicts an exponential increase with time of the perturbation By inspecting equation ()23[) . we see 
that it corresponds to decompactification. 
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the very beginning of the process, and its validity is confirmed afterwards, by proving that 
the chosen ansatz is actually a solution. Of course, the chosen manifold may not provide 
a solution, or the solution may lack some desired properties. 

In looking for a many-fluxes landscape, we have tried the following candidates as com- 
pactification manifold: 

• The complex Riemann curves 

w 2 = P k (z) , 

where w and z are complex coordinates, are the simplest surfaces of known genus 
g > 0. Pk{z) is a fc-degree polynomial with real coefficients, 

P k (z) = z k + a k _ 2 z k ~ 2 H h a x z + a , 

which are the k — 1 deformation moduli of the surface. The coefficients of z k and z k ~ x 
can be removed by rescaling and shifting z respectively. The genus of the surface is 
given by k = 2g + 2 (if k is even) or k = 2g + 1 (if A; is odd). 

These surfaces are not compact, but they can approximate locally a compact surface 
of the same genus. Thus, Einstein equations are to be solved only locally near the 
approximation region. But this ansatz turns out to yield no solutions, not even in 
this approximate fashion. 

• The compact hyperbolic manifolds (CHM) are fundamental domains of nonabelian 
lattices in the Lobachevskian plane, in which the lattices are generated by discrete 
subgroups of SL 2 (M.). By choosing identification of the sides of the fundamental cell 
a compact surface of genus > 2 is obtained, which has constant negative curvature^. 
These surfaces are indeed solutions of Einstein equations, but only if the cosmological 
constant A is negative. This generates a sector of the Einstein-Maxwell landscape 
with no de Sitter states, and therefore will not be considered here. 

The previous examples show that sometimes the appropiate compactification can be elusive, 
maybe because is far more complex than expected, or because it may not exist. So we 
are forced to consider another simple extension of the four dimensional Einstein-Maxwell 
landscape, which is discussed in the following subsections. 

3.1 Multi-sphere compactification 

We consider a J-sphere metric ansatz given by 

j 

ds 2 = e 2 ^ (-dt 2 + dx 2 ) + e 2Mu ' v) (cK 2 + d^ 2 ) . (42) 

8=1 

3 Refs. [371 13H] contain more details on hyperbolic compactifications in cosmology. 
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The metric (142]) represents a manifold of the form 



J spheres 

' A * 1 

K = (A)dS 2 x S 2 x • • • x S 2 = (A)dS 2 x [S 2 ] (43) 

which is nothing but a sector of the Einstein-Maxwell theory in 2 J + 2 dimensions. The 
exponent defines the conformal factor of a two-dimensional cosmological spacetime in 
(t,x) coordinates. The functions ipi depend only on the corresponding coordinates (uj,t>i) 
(but not on (uj, Vj) with j ^ i), and they give a conformal representation of the i-th sphere 
metric. 

The Ricci tensor of the metric (1421) is 



Rxx — —Rtt — 4>tt ~ (frxx i 

RuiUi B, ViVj ^i^Pi ) 

all remaining components being zero. The i-th Laplacian operator is Aj = <9 2 . + <9 2 .. 
The metric ( 142]) should be complemented with the electromagnetic field 

j 

F = ^ Bi(ui, Vi) dui A dvi (45) 

where the magnetic m,, fj- component Bi depends only on the corresponding coordinates 
(ui,Vi) (but not on (uj,Vj) with j ^ i). Thus, the electromagnetic tensor F^ v is analogous 
to (ITU]) , with J 2 x 2 blocks along the diagonal. Maxwell equations ( ITT]) reduce to 

8 Ui y/=g F ViUi = , 3* V=0 F U%V% = • (46) 

Using the volume element prefactor y/—g = e 2< ^ +2 ^ii ^ and the relation F UiVi = e -4 ^^^^ 
we find that a solution of the equations is 

F UiVi = B i = ^-e^, (47) 

where Vi is the volume of the z-th sphere, 

Vi= [ e 2 ^d Ul Adv l} (48) 



and Qi is an integration constant. When integrating the two-form we obtain 

[f = J2w f/' pi du i Adv i = J2Qi = Q, (49) 

so that Q is the total magnetic charge of the configuration, and each constant Qi can 
be interpreted as the magnetic charge contribution of the corresponding magnetic field 
component. 
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The most convenient way of obtaining the field equations is writing the action 

S = — [ ^g(R - 2A - F 2 ) dt dx TT dm dv t (50) 
lo7r / v ' - LJ - 

J 8=1 

in terms of the fields 0, ipi and Bi and then derive the equations from it. The curvature 
scalar of the metric ansatz fl42l) is 



R = 2(<p tt - 4> xx )e-^ - 2 A ^ ■ ( 51 ) 

i=\ 

The Lagrangian density of the electromagnetic field fj45|) is 

j 

F 2 = 2j2B-e~^ z - (52) 
i=i 

Thus, the action specialized to our ansatz is 

S=±-f e^Ek* [(^-^e-^-^A^e-^-A-^Bfe-^] dt d* f[ du, d«, 

i=i i=i j=i 

(53) 

The variation of the action ( 1531) with respect to the vector potential (which determines 
the magnetic field Bi = d Ui A Vi — d Ui A Vi ) gives the Maxwell equations ( 1461) . Varying with 
respect to and ipj gives 



A = ~J2{A^e- 2 ^ + B 2 e-^}, 
i=i 



(54) 

(A, - <#w)e- 2 * -A + A A e- 2 *- + B\ e" 4 * } - e~ 2 * - 2B 2 e" 4 * . 

i=l 

Note that the first equation in (|54p cancels some terms in the second. Now, we recognize 
the Gaussian curvatures of conformally flat metrics with signatures — h and ++; so we 
substitute the constant curvature ansatz implied in the geometry of fC ( 1411]) as we did 
previously ( !T6|) . introducing the constants 

A = (<fo ~ <M e- 2<t> , Ki = -AiV< e~ 2 ^ , (55) 

where A is the curvature of the non-compact part (the cosmological constant of the two- 
dimensional spacetime) and Ki is the curvature of the z-th sphere in the product [S 2 ] J . If 
we finally substitute B t by the solution (jSj) , we obtain 

A-E{*'<-(t) 2 }. (^) 

i=l x / 

X = K 3 -2{9±j. (57) 
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We can substitute equation (1ST)) in (!56|) . obtaining 



(58) 



The relation ( FlTl) between volume and curvature is valid, and transforms equation ( 1571) in 
an algebraic equation for Kf 



47T 



Ki + X — 



(59) 



which has two solutions 



A" 



(±) 



4vr2 

of 



1 ± 



4vr 2 



(60) 



Substituting (160]) in equation (|58|) . we obtain a single equation for A, whose solution can 
be substituted back in equation (|60|) . determining the curvatures. 

It should be noted that if Qj = 0, then equation ( 159)) has a single solution, namely 
Kj = A, which is the limit of the Kj~^ solution when Qj — > 0. 

Now, the usual Dirac quantization condition is 



QjC = 27m j 



with Uj G Z . 



(61) 



which can be justified in the following way. We will use the conformal representation of 
the sphere metric: 



duj + dvj 
Kj cosh 2 Uj 



(62) 



such that the Uj coordinates separate hemispheres (uj > is the northern hemisphere, 
Uj = is the equator, etc.) and Vj are angles mod In. A quantum wavefunction \^ 
defined on the manifold K depends on the coordinates t, x, {uj, Vj}j=i : ...j. The loops 7j in 
which Vj varies along [0, 27r] and the remaining coordinates are fixed can be used to define 
holonomies hj acting on the wavefunction: 



hjty = exp ( te 



AW 



(63) 



where the electromagnetic potential is used as the connection to parallel transport the 
wavefunction values along the loop. It is well known that the potential of a magnetic 
monopole can be defined in two patches on the sphere which overlap at the equator: 



A 



A, 



i=i 




Qi 
Kj\ 
Qj 



- [tanh-Uj — l\dvj if Uj > 0, 
r- [tannic + l] dvj if Uj < 0. 



(64) 
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Here, Vj = as usual, and the magnetic field is 



Bj ~ " KM cosh 2 u 3 ~ V 3 6 A ■ {bb} 

Thus, the potential is discontinuous at the equator, but the discontinuity is given by a 
gauge transformation Xji namely 

•1/ -V <1\,- (66) 

If we move the loop slightly upwards or downwards from the equator, the discontinuity 
in A will leave a different phase on thus violating gauge invariance, unless the phase 
difference is an integer times 2n, that is, equation ( JBTT) . 

Therefore, we can substitute ( ISTT) in the equation satisfied by A (1581) . obtaining 



A = i 
2 



n- 



Note that e 2 is a scale which can be used to measure A and A. In order to simplify the 
formulae, we will assume henceforth that the substitutions -4 — > A and -4 — > A have been 
made. Equation fl67|) takes the form 

L{sMnA X )= A i with (A) = 2 [ JA + E ^2 (i + ^V 1 ~ 2A ^ )] • (68) 

The function L{ s } { n }(A) depends on the signs {s^} of the curvatures involved and on the 
winding numbers {n,j} of the magnetic field. When some rij vanishes, the curvature should 
be taken as Kj = iC ^ = A, the K^ branch is absent. The solutions of equation ( 1681) 
constitute the states of the multi-sphere Einstein-Maxwell landscape. 

The set of integers {rij} will be called a node, while the set of signs {sj} will be called a 
branch of the L^ s .y^ n .-^ function. Fixing a node with all nj nonzero, we have 2 J branches, 
one for each possible choice of the signs {sj}. The number of solutions of the equation is 
different for each A value; in figure [2] it is shown an example with J = 3, where the eight 
branches of the node (m, n 2) n 3 ) = (2,3,1) are shown. The displayed value of A yields 
eight solutions, but is it obvious from figure [2] that the number of solutions vary when A is 
moved upwards, and becomes zero when the {+,+,+} branch is surpassed. The highest 
branch will be always the all— I- branch, and will be called the principal branch. 

Not all solutions of equation f l6"5|) give valid states; for example, curvatures flBD]) can 
become complex if A is too large. The smallest value of A at which some of the pairs of 
curvatures K^ meet is the branching point of the node, and is given by 

Ab = 1 j n • (69) 
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Branches in node (n^ n 2 , n 3 ) = (2,3,1) 




-1.5 -1.0 -0.5 0.0 



Figure 2: The eight branches of the J = 3 multi-sphere Einstein-Maxwell landscape at the 
node (rai,n2,n-3) = (2,3, 1) are shown. Note that each pair of branches meet at a branching 
point, whose horizontal position is the same in all branches, see text. The solid horizontal 
line corresponds to a random value of the cosmological constant A; with this choice, there are 
eight solutions of the equation L(X) = A, yielding a single de Sitter state and seven anti-de 
Sitter states. Note that if A were at 0.5 height, there would be only four solutions, and if it 
were at 1.5 height, there would be no solutions at all. All statistic plots in this paper were 
done using R [39] . 
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This is the maximum positive value that A can achieve when all rij 7^ 0. Curvatures 
can also become negative; this happens when A < in all branches except the principal 
one. Those states are not well defined, but we might define them in detail by replacing 
the corresponding sphere by a CHM. Nevertheless, to keep things simple, this sector of 
the landscape will be deliberately left out, because all its states are AdS. This leaves the 
principal branch as the only source of AdS states with positive curvature in all places of 
the compact part. 

The problem of the stability of the states just found is addressed in subsection 13. 2\ and 
the problem of counting them is the subject of section HI 



3.2 Moduli stabilization 

The next step in the analysis is to determine if the states of the multi-sphere Einstein- 
Maxwell landscape are stable or not. Our approach will follow closely that of subsection [221 
We begin by introducing perturbations £j(t, x) which represent deviations of the curvatures 
found as solutions of the Einstein equations; the perturbed metric ansatz in the Einstein 
frame is 

j 

ds 2 = e 2^)-2E/ =1 ^(^)(_ dt 2 + dx 2j + ^ e 2^-K,^+^(M0( du 2 + ^2) (7Q) 

3=1 

As before, 4>(t,x) represents a 1 + 1 cosmological solution, and ipj( u j, v j) corresponds to 
the metric of the internal spheres in conformal coordinates. The perturbations £j(t,x) 
describe changes in the radii of the internal spheres, and thus they will be referred to as 
the multi-radion fields. 

We will proceed by writing the action (150|) specialized for the metric (ITU]) . The Ricci 
scalar is 



L 3=1 1 

+ 2 fe)«(E&), - afe).(E&). - ^-*-^-^= lS * A .J } 

k^j k^j 



Note that the previous expression involves second derivatives of the multi-radion fields. 
We can replace those terms by first derivatives by integrating by parts in the action. This 
amounts to the following replacement rule: 

te) t / E - & -2fe) t (E&)/ E ^ & > ( 72 ) 

k=i 

and another analogous equation with the x derivatives. 
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The next step is to note that the expression (H7|) for the magnetic field remains un- 
changed, because the perturbations appear in Maxwell equations ( 146|) only as factors de- 
pending on variables (t,x), and thus they can be factored out of the equations. Thus, the 
contribution of the magnetic field to the action is not exactly (1521) but 



2 -4^-4^ 



(73) 



with the same magnetic field Bj = y- e 2 ^ J . 

Finally, we will insert in the action the unperturbed sphere metric ansatz —e 2 ^ j Ajipj = 

Kj. 

Gathering all these ingredients, we obtain the following Lagrangian for the perturbed 
metric: 

L = ^g(R - 2A - F 2 ) 
j 



tt ~ <Vxx ) 



(74) 



3=1 



Now, we can perform the integration in the internal variables {uj,Vj} and thus obtain 
a dimensionally reduced model for the cosmological part together with the multi-radion 
fields: 



S=±[ LdtdxHdu.dv^UllVj) j L^dtdx, 

J .7=1 V J = 1 7 J 



(75) 



with a 1 + 1 Lagrangian 



I 3 = 1 3=1 



(76) 



where the summation in a, e {t, x} is implied, and we have used (see (J28 
{g^) ap dadfi = e 2 ^ dad/3 = e 2<t >{-dt 2 + dx 2 ) , 



R} 2) = -2{<p xx -<p tt )e- 2,t> 



A 
1 



4:71 



(77) 



It is apparent from (17B"j) that all radions are coupled by the dilatonic factors. When 
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exhibiting all fields explicitly we obtain 



L (2) 



2 



(7* 



i=i 



which is the generalization of the Lagrangian previously found for J = 1 , see eq. ( 130]) . Note 
that a substitution rule analogous to (129]) 



Efc), 



3 2 Efc ^fc 



(79) 



has been used in passing from the first line to the second in ( 178]) . The equations of motion 
are as follows; d n ^^- = ^4^- is 



2(E(&«)(E&)/0 + E&)J = -ae 2 * E^*) 



(80) 



and d n &- = ^ are 



w - - 2(ei) a (X)(^)^) + E^w^ 



-fPmti) . (si) 



The absence of perturbations should be a solution of the previous equations. We can verify 
this requirement by substituting £j = in ( 18T]) : 



V al3 <Pape-^ = U'(0) = X. 



(82) 



which is the cosmological solution previously obtained f fT6|) . fl55|) . Using eq. ([77]) . eq. ( 182]) 
is reduced to 



r;,o) =a',-2(^ 



A 



33) 



A Kj X 
1 ~ ~2~ ~ 2 







4) 



that is, exactly equation ( [59]) determining the curvatures. By substituting £ 3 - = in ([80]) 
we obtain 

E^(°) = ° =► ^E 

which is exactly equation ( 158]) . Thus, both equations (183]) . (184"]) are the correct, unperturbed 
ones. 

At this point, we proceed as in subsection 12.21 by fixing the cosmological background 
by means of (182]) . This fixing is equivalent to neglecting the backreaction of the pertur- 
bations £j on the cosmology 0. This being done, <fi is not a degree of freedom anymore, 
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and the variation of the Lagrangian in (1781 with respect to <p is meaningless. This 
background fixing step can also be viewed as solving equation (|80|) to zeroth order. In 
this approximation, the multi-radion field moves in a fixed cosmological background and 
its evolution is determined by equation ( IHTl) . 

We now turn to the analysis of equation (ED). It is a strongly coupled, nonlinear system 
of equations with only one known solution, = 0, which gives rise to the landscape under 
study. There is no hope of finding a nontrivial solution to this system; nevertheless, we 
are only interested in a description of the stability of the trivial solution. 



To this end, we substitute ( 1821) back to equation ( 1811) . obtaining 



A-e- 2 ^[/'(0). (85) 



The derivative couplings appear in a quadratic form. The linear stability analysis allows us 
to approximate the system of equations by Taylor-expanding to first order the right-hand 
side of equation fl85l) and neglecting the quadratic derivative couplings, thus considering 
the much simpler linear system (written in matrix form) 



— e 



-2<j> a/3 



a/3 



(86) 



We have used the symbol £ to denote the J-component column vector of the perturbations 
and the frequency matrix H is given by 



H 



_d_ 



2[{K j -2\)5 jk -X] 



17) 



that is, 



H 



-A 



V 



-A 



-A 
K 2 -3\ 

-A 



-A 
-A 



Kj-3\) 



The linear stability condition is therefore that all eigenvalues of the matrix H should be 
positive. We will call k the minimum eigenvalue of H, so that the stability criterion is 
simply 

k > . (89) 

In order to complete the linear stability analysis we should justify the neglecting of the 
derivative couplings. These non-linear terms have a generically negative sign9, as opposed 
to a positive sign characteristic of a dissipative force, and thus they could be interpreted 
as an "anti-dissipative" force. It is legitimate to ask if these non-linear terms can spoil the 
linear stability of the solution. 



4 Note that the derivative couplings appear in a quadratic form having all negative eigenvalues except 
for one, see the paragraph following equation (|154l) in appendix [X] 
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It turns out that these nonlinear terms do not spoil the linear stability criterion ( )89|) 
as long as the amplitude of the perturbation is sufficiently small. The magnitude of the 
threshold and the corresponding heuristic argument leading to these conclusions, not being 
central to this discussion, have been placed in appendix |A] 

The exact computation of the spectrum of H is not possible in the general case. There- 
fore, we cannot derive a formula k(A) to quickly establish the stability of a state. As a 
result, the computation of k should be done numerically on each individual state. Never- 
theless, we can obtain some general stability results by computing the determinant of H, 
a task which can be done as follows. 

Firstly we note that if all Kj are equal, the determinant of H would be (we drop from 
now on the unimportant factor 2 of H) the characteristic polynomial of the matrix 



B 



/3A A ■•• A\ 
A 3A ■•• A 

^A A •■■ 2,\) 



(90) 



This matrix has an eigenvalue 2A with degeneracy (J — 1) and a simple eigenvalue (J + 2)A. 
Its characteristic polynomial is 



j-i 



det (ATI - B) = [K - ( J + 2) A] 2A) . (91) 

3=1 

Specializing the variable K at a single curvature value Kj would give the determinant of H 
if all curvatures were equal to Kj. This is not the value of the determinant we are seeking; 
but we can form a permutation-invariant superposition of all those expressions: 

.. j j 
det H = - [Ki - (J + 2) A] \\{Kj - 2A) . (92) 
i=i j=i 

3+i- 

The factor 1/J comes from a normalization condition. This turns out to be the correct 
expression for the determinant of H, and it is straightforwardly transformed in the char- 
acteristic polynomial of H . 

1 J J 
det(# - ill) = -^[Jfj - fi - (J + 2) A] Y[{Kj -fi-2X). (93) 
i=i j=i 

Nevertheless, the computation of its roots is not possible in general. 

Based on expression ([931 , it follows that whenever A < the determinant det H is 
positive, and furthermore det(H — //I) cannot vanish at a negative value of fi, and thus the 
stability eigenvalue should be positive. As a result, all AdS states of the model are stable. 
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Another case worth investigating is those states which have at least a vanishing quantum 
number rij = 0. The corresponding curvature is not defined in this case, because 

equation ( 159]) is linear and it has only the solution = A. Thus, we can substitute 

K j = A + 5j (5j > 0) in det H. Assuming we can vary independently A and 5j, we can 
expand det H for small A values: 

j j 

det#^> \\5 l -2\^\\5 J + 0{\ 2 ). (94) 

i=l i=l j=l 

The previous expression shows again that a negative value value of A cannot make this 
determinant to change sign. A positive value indeed can change the sign in the determinant, 
and this indicates that de Sitter states are likely to be unstable. Of course no general 
statement of this sort can be formulated, because this depends on the magnitude of A as 
well as on all 5f for sufficiently small A > and fixed 5j, the determinant can be positive. 
But if a single 5k = 0, then the determinant reduces to 

j 

detH -2AJJ5, + 0(A 2 ), (95) 

3=1 

which is certainly negative for A > 0! Thus, we conclude that an odd number of eigenvalues 
of H have changed their signs and the state is unstable. We can suspect that in this case 
a single eigenvalue has reversed sign, because if we would take two vanishing 5k the sign of 
the determinant would again be positive. We can see that if all 5j = 0, H has a completely 
negative spectrum. Thus we can expect that the vanishing of each 5j changes sign of an 
eigenvalue, and thus all states with some rij = are unstable. 

The previous reasoning is heuristic, because we cannot assume that A and 5j vary 
independently. They depend on the discrete numbers {m, • • ■ , nj} and thus the stability 
criterion should be validated numerically. Nevertheless, heuristics works in this case. As 
it is shown in the following subsection, all states with a single rij = are unstable, at least 
in all searches we have carried out. 

We will close this subsection by giving a perturbative argument showing that all low-A 
states lying in all non-principal branches are unstable. This way, the principal branch 
remains as the only source of AdS and stable dS states. We begin by splitting the H 
matrix ([58]) as follows: 

H = 2 diag^ -2\,-..,Kj- 2A} i=1) ... , j - 2XU , (96) 

where U is a J x J matrix filled with ones. If A is small, we can consider the diagonal matrix 
as a "dominant" term and the off-diagonal terms as a perturbation. The perturbation is 
permutation-invariant, and thus all eigenvalues of the diagonal matrix are shifted the same 
amount. The perturbative stability eigenvalue is 

k = 2 [min Kj - 3 A] . (97) 
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If the minimum curvature is taken from the (+) branch, then Ky~ — > \ and the state 

has a chance of being stable. But if some curvature is taken from a (— ) branch, then we 

have Kj ^ A ^°> A and the corresponding eigenvalue is negative. But then k should be 
negative also, showing that no matter how small A might be, if the state comes from a 
non-principal branch, it will be unstable. 

The previous argument does not rule out the existence of higher A stable dS states in 
non-principal branches, but our numerical searches have not found them. 



3.3 State searching in concrete examples 

Once the analysis of the model is reasonably complete, we should ask for concrete exam- 
ples where we can exhibit some states and their associated magnitudes such as reduced 
cosmological constant A, curvatures K t and the stability eigenvalue k. 

We repeat here the relevant equations for the reader's convenience. Given a J-tuple of 
integers {rti, • ■ ■ , rij} we compute the solutions of equation fl68l) . which is 



J {sj},{nj} 



(A)=A, with L {s . }){n . } (X) 



1 m 

3=1 J 



2Xn) 



(9£ 



All solutions for each branch {si, • • • , sj}, where the Sj are signs ±1, should be computed. 
The corresponding solutions must be real, and all its curvatures must be real and positive: 



K 



71: 



2Anf 



(99) 



If a single curvature turns out to be real and negative or complex, then the state should 
be discarded. If all curvatures are positive, we form the frequency matrix H (dropping 
the unimportant factor 2 which appears in ( 188]) ) and compute its minimum eigenvalue k, 
called the stability eigenvalue of the state: 



H 



{Kx-3\ -A 
-A K 2 - 3A 



-A 
-A 



\ 



V 



-A 



-A 



K, 



3XJ 



The stability condition is simply k > 0. 

Thus, the searching method has the following steps: 

1. Choose J and A, the parameters of the model. 

2. Choose a set of integers {n\, ■ ■ ■ , rij}. 

3. Solve equation (1981) for A. 



mm /i . 

/^Gspec(_ff) 



(100) 
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4. Compute the curvatures (|99|) associated with the solutions found and accept the state 
if all curvatures are real and positive. 



5. Compute the stability eigenvalue and flag the state as stable or unstable. 

6. Go to step 2 until some bounding search criterion is met. 

The choosing of the integers can be done in various ways. In the simplest models with 
we can scan a large square in the (ni, 712) plane in a brute-force search. In this way 
we cannot miss any state. The symmetries n% n 2 and rij — > —rij allow us to restrict to 
n 2 > n i > 0. This brute-force method cannot be used for higher values of J. In those cases, 
we should generate states randomly in an efficient manner; but prior to the discussion on 
how this is done, we present some results in the J = 2 case. 

We have chosen two values of the 4D cosmological constant A, 0.01 and 0.005, in order 
to exhibit how the lowering of A causes the proliferation of states. In figure [3] we plot a 
point in each place of the (7x1,722) plane where a state has been found; of course, we have 
four branches for searching, so some states overlap here. States with negative and positive 
2D cosmological constant have been separated, so that we can see AdS states at left panels 
and dS states at right panels. Stable states have been drawn using circles, and diamonds 
for unstable states. Note that all AdS states are stable, while most dS states are unstable. 
The colors are related with the magnitude of the cosmological constant as shown in the 
legend of each graphic. 

Features of these models which can be seen in figure [3] include: 

• All AdS states are stable. 

• There is no AdS states with n\ = or n 2 = because in those cases there is no 
curvature associated with this branch. 

• All dS states with rii = or n 2 = are unstable, as expected by the heuristic 
argument exhibited at the end of the previous subsection. 

• There is a curve which limits the existence both of dS and AdS states. The form of 
this curve is easily computed by substituting A = in equation ( 1981) with all positive 



Minkowski states, which ideally may be present, lie on this hypersurface which will 
be called branching hypersurface (curve in the J = 2 case). In practice, Minkowski 
states require fine-tuning of A and thus they are generically absent. 

5 The two-sphere Einstein-Maxwell lansdcape is a sector of the six-dimensional Einstein-Maxwell not 
considered in Refs. [?T1 155] . 




(101) 
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• All dS states are located above AdS ones, and they are located near the branching 
curve. The discrete states are located on a multi-branch surface whose branches meet 
at the branching curve, hence its name. Thus, states of A near zero from either side 
are located near this curve, which would contain, if present, the A = states. 

• All stable dS states are located near the branching curve, and only there, but closeness 
is not enough for a state to be stable, as will be seen below. 

A feature of these models which cannot be seen in figure [3] is that all AdS states come 
from the {+, +} branch (which will be called the principal branch). This is so because the 
solutions in the remaining branches have the iC ^ curvatures negative, as can be seen in 
formula ( 1991) . Thus, only the principal branch is a source of AdS states. 

In contrast, dS states can come from each of the four branches, but each branch brings 
in states with very different properties. For example, all dS states near the branching curve 
come from the principal branch. In particular, all stable dS states come from this branch. 
The bunch of dS states lying in the reddish square in the right panels of figure [3] come from 
the {— , — } branch, and they are "most" unstable in the sense that they have the greater 
(in absolute value) negative stability eigenvalue. All these features can be seen in the left 
panels in figure HI These left panels show cosmological constant versus curvature, and they 
exhibit clearly the different nature of the states coming from different branches: 

• The states coming from the principal branch (bullets) form the core of the figure, and 
at the bottom of this figure are located the stable states, all of them coming from 
the principal branch. Note that stability does not mean lower cosmological constant! 
In the right panels of figure 0] we can see that stable states are mixed with unstable 
ones in the cosmological constant value distribution. Nevertheless, all lowest-lying 
states are stable. 

• The states coming from the {+, — } and { — ,+} branches are at both sides of the 
triangle-shaped graphic shown at the left panels of figure |U All of them are unstable, 
but the values of the cosmological constant in this subset range from lowest to highest. 

• The states coming from the {— , — } branch are located at the cusp of the triangle 
(figure HI left panels) and all of them have the highest values of the cosmological 
constant and also the highest (in absolute value) negative stability exponents (see 
figure [5]). In particular, the n\ = 0,n2 = state, which is the cusp of the triangle, 
belongs to the { — , — } branch. 

In the right panels of figure H] we can see the cosmological constant distribution of dS 
states. The contribution of the branches can be seen as different peaks; while the principal 
branch contributes with the stable states and other unstable distributed in the lowest 
range, the {+, — } and { — , +} show two peaks in the low and middle range, and the bulk 
of the { — , — } states are relegated to the high range. 

Finally, in figure is shown the behaviour of the cosmological constant versus the 
stability eigenvalue. Its almost linear relation can be seen to be dispersed in branches, 
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J = 2 Landscape (AdS states, A = 0.01) 



J = 2 Landscape (dS states, A = 0.01 ) 




J = 2 Landscape (AdS states, A = 0.005) 



" = -0.03372696 



J = 2 Landscape (dS states, A = 0.005) 



! 
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♦ 



♦♦♦ 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 



Figure 3: Landscape examples with A = 0.01 (top panels) and A = 0.005 (bottom panels). 
AdS (left) and dS (right) states are shown, using a circle for a stable state and a diamond 
for a unstable one. A state is drawn if a solution has been found for A with real and positive 
curvatures. The magnitude of A is reflected in the color of each symbol. Finally, all dS states 
lie above AdS ones, and meet at the branching curve, beyond which the landscape has no 
states. 
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Figure 4: Landscape examples with A = 0.01 (top panels) and A = 0.005 (bottom panels). 
Left panels show the variation of the cosmological constant A with curvature. These triangle- 
shaped graphics show structures which the states seem to follow, and clearly separates the 
different branches by the A values they provide. Right panels show the A distributions of dS 
states, whose peaks come from the different branches. The only branch that has not a peak 
associated to it is the principal branch, which provides all stable dS states. 
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Cosmological constant vs stability eigenvalue (J = 2, A = 0.01 ) 



Cosmological constant vs stability eigenvalue (J = 2, A = 0.005) 




V 
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Figure 5: Landscape examples with A = 0.01 (left panel) and A = 0.005 (right panel). 
Dependence of the cosmological constant versus the stability eigenvalue is shown, with an 
almost linear behaviour. Stable states are located to the right of the vertical dashed line, 
which is thus the onset of stability. 



which are the same structures showing up in left panels of figure HI Here, we can see that 
the states with lower cosmological constant have a "less negative" stability eigenvalue than 
the states with higher A, which include the {— , — } states, as said above. 

Thus, the different branches provide very different states. Among them, the most 
interesting seem to be those coming from the principal branch, because they include AdS 
as well as stable and unstable dS states, which are the ingredients we need to construct a 
toy model of a multiverse. 

For J > 2, we can exploit the fact that stable dS states are near the branching hy- 
persurface (it is a curve only for J = 2) and design a sampling method which looks for 
states in (n 1; • • • , nj) space whose Voronoi cell (which is the cube with its center at the 
point in question) intersects the branching hypersurface. Those states are called secant 
states [Ull E], and all states near the branching surface belong to this category (but the 
reciprocal is false, that is, a secant state may not be near the branching surface if J is 
large enough!). Thus, we can sample the principal branch by sampling the secant states. 
In this way, we always find a state with a fair chance of being a true state of the model 
in the principal branch, that is, we have an efficient sampling method, much better than 
brute-force node enumeration or completely random node sampling. 

The sampling of secant states is simple: we choose a uniformly random direction in 
J-space and find the point of intersection between the ray having the chosen direction 
starting from the origin and the branching hypersurface. This intersection point belongs 
to the Voronoi cell of a single secant state, which we find by rounding the coordinates of 
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the intersection point. Once we have the state, we solve the equation for A in the principal 
branch and follow the steps detailed above. 

In the following section we use this sampling method to obtain a sample of the cosmo- 
logical constant distribution which can be compared with an approximate formula to be 
obtained below. 

4 State counting 

As we have seen in the previous section, the states of the multi-sphere Einstein-Maxwell 
landscape can come from different branches, and the richest of those branches is the prin- 
cipal one. In this section we turn to the problem of counting states on this branch. Our 
main aim is to compute the distribution of cosmological constant values in this branch. 

As discussed above, only the principal branch can have both AdS and dS states, both 
stable and unstable, and moreover the large-J sampling method is especially adapted to 
the principal branch. Therefore, we can obtain samples to compare with the approximate 
formula to be obtained in subsection 14.21 below. 

4.1 Counting states in the principal branch 

Given a node n = (ni, • • ■ ,nj) in J-space, there exists a state in the principal branch of 
the multi-sphere Einstein-Maxwell landscape with cosmological constant A if the equation 



has a solution in A. The principal branch L n function is given in (198]) with all positive 
signs Sj — +. In the principal branch there are no states with some rij — 0, and thus the 
function L n (X) has a maximum at A = and it is monotonically decreasing from A = 
to A = Ab, where At,, see equation flo^]) . is the branching point of the node n. A dS state 
exists therefore if A is between the two extremal values of the L n (X) function: 



Equation ( I103p is the existence condition for a dS state at node n in the principal branch. 
The corresponding equalities define two surfaces in node space: The branching surface 



L n {\) = A 



(102) 




(103) 




(104) 



whose integer points, if any, have vanishing A, and the limiting surface 



-^n(Ab) 
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(105) 
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which signals the end of the principal branch. In terms of a characteristic function, the 
existence condition is 



J 1 if a state exists at n, 
X [LMLnm (A) = | q . f a gtate doeg ^ ex . st ^ n (106) 

The previous existence condition should be supplemented with the stability condition ( 189|) . 
that is, K(n) > 0, where K(n) is the minimum eigenvalue of the stability matrix H at node 
n if a state exists there. This stability condition can also be ascribed to a surface (the sta- 
bility surface) signaling the stability threshold. The stability surface should be comprised 
between the branching and limiting surfaces; unfortunately, the analytic expression of it 
cannot be found for general J > 1. As a consequence, we will represent this condition by 
adding the factor 6(k(ti)) to the existence condition. Therefore, the exact number of stable 
states with given A in the principal branch of the multi-sphere EM landscape is 

Afj(A) = X[i„(A b )A(o)](A)e(K(n)). (107) 

n&Z J 

The exact evaluation of the previous expression is possible only for J = 1 as is showed 
in subsection 12. 2\ equation (flTj) . We repeat it here for the reader's convenience, omitting 



M(A) 
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We will shortly turn into the approximate evaluation of ( 110711 . But prior to that, we need 
to grasp some general ideas on the structure of the stable state set which we are willing to 
count. 

We will emphasize two main aspects: asymptotic hyperplanes and state chains. 

Asymptotic hyperplanes Both equations (110411 and (I105f) corresponding to the branch- 
ing and limiting surfaces, and likewise the stability surface, have asymptotic hyper- 
planes located at 

\rij\ = —= (for fixed j), \nk^j\ — > oo , (109) 



X A 

and thus all states are restricted to the region 

1 



nj \> — = v (l<j<J). (110) 



Therefore, all dS states should have a charge greater than u as a necessary condition. 

Note that no integer rij can be equal to u while preserving the existence condition 
unless all the remaining integers n^j are infinite, hence the name "asymptotic". It 
is easy to see that the corresponding states would have A = and all curvatures 
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vanishing except for one, and thus they would not represent compactified states. 
Moreover, the associated stability matrix has J — 1 zero eigenvalues (see eq. (HOOD ), 
and thus these states are only marginally stable. These properties suggest that they 
should be excluded from the landscape. 

State chains There is also a natural upper bound on the charge, which can be obtained 
by considering the following straight line in flux space: 



nj-i 



J 



A 



n,j G IR (free parameter) . 



'Ill) 



The previous line is asymptotic to the branching surface, in the sense that it satisfies 
equation (11041) when rij — > oo. This line do not contain nodes because the quotient 



~~ic~~ is generically not an integer. But we can slightly modify the previous line: 



nj-i 



J-1 



A 



Vi , rij G K (free parameter) 
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This modified line can contain valid states if rij lies between \vq\ and z^, where 
is the intersection height with the branching surface: 



V2 



^Av\ -(J-1) 



113) 



All states on the line above z/ 2 are beyond the branching surface. Therefore, we have 
the bound Vq < rij < v-i (1 < j < J). All states which might happen to lie on this 
line are said to form a state chain. 



It should be noted that when 



j-i 



coincides with the integer u±, then z/ 2 diverges, 



which at first sight would be interpreted as an infinite dS state chain of ever decreasing 
A. But those nodes in the chains have no states, as can be seen by explicitly writing 
equation (1 102[) for the nodes in the line f 1 1 1 2 p . and look for solutions with small A 
and large rij\ 



2A = JA + ( J - 1 



1 + vl - 2\v{ 1 + y/1 



2\n 2 j 



v 



We can consider A -\, but A 
equation (I114p can be rewritten as 

2(J- 1 



i 

'In 2 , 



(114) 

i n j 
is not true, because nj is large. Therefore, 
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1 - 2\n 2 j 



n. 
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The right hand side of (11151) is positive, thus a solution to equation f 1 11 5 j) can exist 



only if A 



j-i 



is strictly positive, which leads to 



J-1 



A 



> 



J 



A 



(116) 
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that is, if y * s an i n ^ e g er then there is no solution to equation f lll5p . As a 
consequence, z/2 can be made as large as we want by fine-tuning A but it is never 
infinite. This argument shows that all state chains are finite. 

This discussion on the state chains can be generalized to other asymptotic affine 
manifolds that the branching surface can have. For instance, asymptotic hyper- 
planes f )109p . as we have seen above, are likewise devoid of states, but there are close 
hyperplanes (having nodes) each one containing a replica of a (J — l)-dimensional 
landscape. 

We will now consider the evaluation of the number of stable states A/j(A). Firstly, we 
get rid of the sign degeneracy 2 J , which is always trivially present. Secondly, we invoke 
the permutation symmetry, which allows us to arrange the integers rij in decreasing order. 
The corresponding permutation degeneracy is J! except on those nodes having repeated 
components. This difference will be ignored for simplicity; we will see below that it will be 
of little importance in the small-A region. Thirdly, we consider a node n = {rii, ■ ■ ■ ,nj} 
with ni > n 2 > ■ ■ ■ > rij and the corresponding equation for the existence of a state lll02p : 



L n {\) 



JX 



j= i 



A. 



(117) 



The L n (X) curve has a branching point Ab given by fl69l) . that is, Ab = Any approxi- 
mation method we might wish to apply on the L n (X) curve should respect this branching 
point in order to accurately represent the existence condition. In particular, we cannot 
assume A <C Ab- But we do have rii > rij = miiij{nj}, and in the case rij <C ni, we can 



assume A i and write 

2n 2 j 



which leaves equation 1 1 11 7ft as 



Kj(X) 



;ii8) 



j-i 

(J-l)A + X>i(A) 

j= i 



A- — . 



;ii9) 



Equation 1I119P represents the solutions of a landscape in which the J th curvature has been 
removed, and the cosmo logical constant A has been replaced with A — \. We can now let 



nj run from |_^J + 1 through the diagonal node having rij — [ 
recurrence law 

AO (A) w J\ ^-A A 

The previous formula is valid under the following conditions: 



J , thus obtaining the 



(120) 
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The fraction of states with repeated components is small. 



• The cosmological constant A of the states included is small, so that equation ( 11 18ft 
can be valid. 

The states near the asymptotic hyperplanes will satisfy the previous conditions more accu- 
rately, so that the first terms in the sum ( I120p will be more precise than the terms near the 
diagonal. The latter states will fail to satisfy the strong inequality n; < ni. This means 
that the low-lying (that is, small-A) states will be taken into account, but the formula can 
miss or overcount some high-lying (high-A) states. 

Equation f ll08j) triggers the recurrence relation, the first consequence being 



A/2 (A) 



LvlJ 
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m 2 



2a/2 



■3 J A 



1 

rn 2 



(121) 



The previous equation can be approximated by a smoother version by removing the floor- 
ceiling functions inside the sum: 



JV 2 (A)«2(l 
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The previous formula can be refined by cutting off the chains whose lower ends pass the 
diagonal and counting the states on the diagonal accurately. This "diagonal corrected" 
formula is to be used in figures [6] and [7] below. 

Further simplification can be achieved by isolating the first term (which carries the 
discontinuities) and estimating the remaining sum by means of an integral: 



JV 2 (A)»2 1 




dx 



A 



(123) 



Formulae ( 11221) and f l 1 231) show clearly the effect of state chains as discontinuities at integer 
values of -7=. When -7= approaches an integer from below, a very long state chain develops 
which increases dramatically the number of states. When A is reduced, the number of 
"bulk" states, that is, those not in the chains, increases as reflected by the well-behaved 
integral contribution, which for A 1 is 



dx 



A<1 



V2 



\/A 



A 



A A 



(124) 



Figure provides a good example of the performing of equations ( 1121 j) . (1 122j) and 
( I123p . These formulae are to be compared with brute-force determination of the number 
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Number of stable states versus A (J=2) 



Number of stable states versus A (J=2) 
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Figure 6: Number of stable dS states as a function of A for the two-sphere Einstein-Maxwell 
landscape. Left panel shows the strongly discontinuous variation of the state number when 
^ crosses several integer values. Bullets are brute-force computed state numbers, thin 
green line is the outcome of formula (I121D . thick solid line is formula (I122|) (with diagonal 
corrections) and thick dashed line is formula (| 123|) . Right panel shows an amplification of the 
small rectangle shown in left panel. Data are to be interpreted as before, with the addition 
of hollow bullets, which mark the data displayed in left panel. Simplified formulae seem to 
have an averaging effect on the lattice details, which are reflected in the fast-varying nature 
of the discrete formula and brute-force data. 
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of stable states in the corresponding models. The discreteness of the lattice induces strong 
fluctuations in the actual number of states, which is well represented by formula f !123j) . 
provided we interpret it as an average behaviour. 

While figure |6] emphasizes the strongly discontinuous nature of the state number, we 
can also show the steady increase in the state number by avoiding the discontinuities. This 
can be done, for example, by sampling landscape models with half-integer values of 
These samples never encounter large state chains and thus a regular, well-behaved curve 
emerges, very well described by the formulae just obtained. This smooth component of 
the state number is illustrated in figure CD 




Figure 7: Number of stable dS states as a function of A for the two-sphere Einstein-Maxwell 
landscape. Only half-integer values of have been considered in this plot, in order to avoid 
the spikes shown in figure EJ Bullets are brute- force computed state numbers, thin green line 
is the outcome of formula (|12ip . thick solid line is formula (|122j) (with diagonal corrections) 
and thick dashed line is formula (I123p . A smooth behaviour is observed, showing a very 
good agreement between approximate formulae and numerical searches. The narrow vertical 
rectangle is at the same position as the rectangle shown in figure [6] (left panel). 

We will close this subsection by summarizing the properties of formula H123[) as follows: 
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• It accurately captures the spikes in the state number when ^= approaches integer 
values from below. These spikes come from the presence of very long state chains in 
this regime. 

• It correctly represents the main behaviour of the state number in a generic sense, 
that is, when ^ is not near integer values. 

• We can interpret formula (11 23ft as an average behaviour which turn the fine details of 
the lattice into a smooth profile while taking into account the main discontinuities. 

• Finally, the approximation formulae seem to be missing some states. The reason 
for this is the approximation we are using to count stable states: Equation f 1 1 2 j) 
implies using the J — 1 stability criterion to count stable states in the J model, 
which introduces the error. In the following subsection we will see that the missing 
states are located near the discontinuities of the density of states, and thus they 
correspond to relatively high values of A. This is precisely the condition which makes 
(11191) to break down, so this behaviour was to be expected. 



4.2 Small cosmological constant distribution 

In this subsection, we will call A(n) the 1+1 cosmological constant of a stable dS state at 
node n (assuming the state exists), and we will denote by p a fixed value to be compared 
with A(n). With this in mind, we define the distribution function of p in a given multi- 
sphere EM landscape as the number of stable dS states whose A(n) value does not exceed 
P- 

Vj( P , A)=J2 xiL MLnm W e(K(n)) e( P - ah) . (125) 



The derivative of Qj(p,A) with respect to p is the density of states of the model 

Wj( ,,a) = w^. (126) 

As a result of the discreteness of the landscape, A values are drawn from a discrete set, 
and thus ttj(p, A) is a stepwise- varying non-decreasing function of p, while uj(p,A) has 
Dirac deltas at the values of p coincident with actual A(n) of existing states at n. The 
amplitudes of the Dirac peaks are given by the degeneracies of the corresponding states. In 
this subsection we will obtain some analytic approximations of the density of states in the 
regime of small A, and we will use the expressions thus obtained to study the A spectrum. 

We will denote by A max the maximum A value a state can have. Clearly, if p > A max 
then 

fi J (p,A)=A0(A). (127) 

Analogously, we denote by A min the minimum A value a state can have. It is also clear 
that, if p < A min , then 

fij(p,A)=0. (128) 
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Thus, the interval [A mm , A max ] is the support of the density uj. 

The upper bound A max can be computed as follows. Let us consider the gradient of 
the function A(n) computed as if the components of n were continuous variables. We can 
derive equation (II 17ft implicitly: 



Jd n .\ + J2{^jdn 3 K t + d nj Xd x K t \ = . 
i=i > 



(129) 



It follows 



0..X- d "'?> =1 ' ^ . (130) 

The denominator of (I130p is clearly negative, and thus the gradient is always pointing 
from the branching surface to the limiting surface. Right at the branching surface A = 
and the gradient is always infinite. Right at the limiting surface all components of the 
gradient vanish except that of maximum tij. Thus, the gradient is always pointing towards 
the diagonal, except right at the diagonal, where it points towards the origin. From this 
gradient configuration we conclude that the maximum value of A is achieved at the cusp of 
the limiting surface. But this point is not in the stability window, and thus the maximum 
A should be achieved at the onset of stability along the diagonal. But right on the diagonal 
the stability matrix H is permutation-invariant and we can compute exactly its stability 
eigenvalue, which is 

K A , Ag = K-(J + 2)\, (131) 

where K is the common value of all curvatures on the diagonal of flux space. But then, 
equation 1 1 11 7ft reads 

\ = l[\ + K\, (132) 
which allows us to eliminate K and gives the exact diagonal stability condition: 

x< WTrr x — (133) 

Thus, we have exactly computed the maximum A value. 

Things are far more difficult when we address A mm . We know that the minimum will 
be close to the branching surface, but its exact position is unpredictable in general. The 
J = 1 case is easier, because the landscape is a single state chain. In this case, equation 
(11321) is exact, so we can obtain the dS spectrum as 





a( " ) = 2 (n-v^- a JV^- a ' (134) 

The end state of this chain is at node n max = LttjJ > an d this implies 



1 





'7VL5F- A - (135) 
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In this case, it is possible to give an exact answer to the minimum A value. When ^= is 
an integer, then A m ; n = 0, but, as we know, this is not the generic situation. 

Generalizing the result (11351) is difficult. We can argue as in the previous subsection 
and approximate the J = 2 case by the J = 1 case just considered. Then we can assume 
that the longest state chains will host the minimum A states at its end nodes. These end 
nodes have an approximate J = 1 spectrum which can be computed by replacing A in 
fjT32]) (with J = 1) by A - ([^J + l)~ • Now the same substitution can be performed 
in formula (I135p . giving a fairly cumbersome expression of nested fractions, square roots 
and floor functions, and hence it will be omitted. This formula will approximately give the 
minimum A provided the corresponding state is located at the end of the longest chains. 
The minimum-A state can also be located among bulk states, but this is non-generic, as 
we will see shortly, because long state chains are generically low-lying states. So such a 
formula can be trusted, but it is neither exact, nor a bound, but rather it is an approximate 
expression for the particular (if generic) case when the min-A state is located at the end of 
the longest state chain. 

Similar arguments can give analogous (but much more complex) expressions for higher 
J, with the same caveats as before. 

Nevertheless, A min is not quite relevant for our purposes, because a precise computation 
of it requires taking into account even the finest details of the lattice, and thus no continuous 
approximation can yield this value. Instead, we will be interested in an approximation of 
the density of states, whose expression will allow us to estimate A m i n in more familiar terms. 

Let us consider first the case J = 1. By inverting the relation X(n) = p, we obtain 



n p = , (136) 



so that the condition X(n) < p can be rephrased as n > n p . Thus, the distribution function 
Qi(p,A) is simply the number of integers between n p and n max , that is, 



Q 1 (p,A) 
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2A-p 



(137) 



supplemented with conditions H128[) and (I127p . This exact result can be approximated by 
a continuous function by simply omitting the ceiling functions. Doing this and taking the 
p-derivative afterwards, we obtain 

(2A - p) VA - p 

Equation (I138p is the J = 1 density of stable dS states. It has the following properties: 
M(A) w / u 1 (p,A)dp=(l-^^ ' 



'A 

11 r- (139) 

If 12 + 12 tan" 1 A= - 3tt - 7^2 

<A> - = jm i p " l(p ' A) dp = A ■ 
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That is, it is consistent with equation (jl08p . and the mean value of the density is around 
0.7167|. This curve has a jump discontinuity at p = ~, which is the upper limit of its 

support, and the position of its maximum, which is ^f. Figure M illustrates this density 

compared with the actual spectrum of a J = 1 model. 



Density of states and X spectrum (J = 1 , A = 1 ) 
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X/A 
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Figure 8: Density of states of the J = 1 landscape (thick line) compared with the actual A- 
spectrum for A = 10 -8 . This model has 572 stable states, and the density shown is normalized 
to unity. The agreement between discrete data and continuous density is complete, in the 
sense that the histogram, which approximates a continuous curve when the spacing between 
neighboring states is much smaller than the bin width, accurately fits the approximation 
uJi(X, A). The construction of the estimate A m ; n is also shown: the small triangle located 
at the origin has area 1/572 in this model, and its vertical side marks the position of the 
minimum-A estimate. This value of A allows for a A = state, which is isolated by the bins 
used in the histogram. 

In figure [HI is also displayed a naive estimate of A m i n , defined as the abscissa A m i n which 
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encloses area 1 under the density's graph: 

/ co 1 (p,A)dp=l. (140) 
J o 

Assuming that A min is small enough, we can approximate Ui(p,A) = + 0(p 2 ) and 
obtain 

A min = 2V2 At . (141) 

Figure |9] shows this estimate versus the exact minimum. Of course, this estimate does not 
provide the true minimum: it has neither zeros nor peaks, but it grows at the same average 
rate. 




Figure 9: Comparison between the exactly computed minimum-A (equation (|135h ) and its 
estimate A m i n (equation ()14ip ). The latter seems to be an upper bound of the former, closely 
following the decreasing of the worst-case minimum as A decreases. The spikes shown by 
log A m i n are the values for which -j= is an integer. This continuous-density-based estimate is 
not accurate because the value of A m i n is dictated by the finest details of the lattice and not 
by the continuous density wi(A, A). 
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We will now consider the J = 2 case. The recurrence relation (11211) extends to distri- 
butions and densities as well, and thus we have the approximation 



w 2 (p, A) 



LvlJ 
E 



u x [p, A =■ 



(142) 



We interpret this equation as a decomposition of the J = 2 landscape in a superposition of 
several J = 1 landscapes, which are state chains. Of course, the previous equation will be 

only valid for small p. The maximum value of m in the sum dH2D, which is [,/f J, gives a 



maximum effective cosmological constant A e fj jmax = A — I 
the corresponding state chain is 



A J 



Aeff.max A 

2 4 ' 



and the maximum A of 



(143) 



in contrast with f |133j) . which gives A max = jr. Thus, we see that this approximation gives a 
wrong maximum A value. The origin of this discrepancy is the stability condition, because 
the superposition of state chains extends the validity of the J = 1 stability criterion to 
J = 2, and this is true only for small A. 

An example of this density of states compared with actual A-spectrum data is given in 
figure [10J We can see that the histogram shows a peak near the origin, and the density of 
states extends its support to j instead of the correct j value. In the logarithmic version of 
the histogram, we can see the first peaks resolved enough, and the correctness of the state 
chain approximation in the low-A region. Only a few peaks get resolved; the remaining 
peaks merge in a bulk distribution whose approximation computed from ( I142j) is incorrect 
in the high-A region. We will not need this bulk distribution here. 

The mean value of the u>2 distribution can be directly computed from equations (1 1 4 2 [) 
and ( I139p . We denote the summation interval as /(A): 



and then we have 
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Cosmological constant (stable, J = 2, A = 1 e-04) Cosmological constant (stable, J = 2, A = 1 e-04) 
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Figure 10: Comparison between a brute-force-computed A-spectrum (histograms) and the 
approximated density of states (thick line) in a J = 2 model of the Einstein-Maxwell land- 
scape (A = 10 -4 ). Left panel: The ordinary histogram shows a narrow peak near the origin 
(most of them come from long state chains) and a tail of bulk states. The jagged density of 
states accurately accounts for the first few peaks, but it fails to describe the high-A values. 
Right panel: The logarithmic histogram shows the resolved structure of the first peaks, well 
described by the density of states. This histogram also accumulates the bulk states in a 
single, broad peak. The corresponding bulk part of the density of states is shifted to the 
right. 
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In (I145p . the constant £ is the prefactor of A in the formula for (A)^ 1 appearing in equation 
(11391) . An illustration of the general behavior of the mean cosmological constant is given in 
figure CLU The average A diminishes towards zero when ^= approaches integer values from 
below, as a consequence of the development of large state chains. Formula f !144j) is not 
diagonal-corrected (see above), which gives small unevenly-spaced jumps. It is compared 
with brute-force-computed averages, which fluctuate because of lattice details. The global 
decreasing of the mean value as A decreases when ^= is half-integer is shown also in figure 
im (right panel) . 

Figure [TT] also shows the minimum A, computed using formula (11411) with the longest 
chain of J = 2 models, compared with brute-force-computed minimum values. The fluc- 
tuation here is caused by the unpredictable nature of the minimum, which can be located 
at any point near the branching curve. When ^= approaches an integer, the very long 
state chains are mainly formed out of low-lying states, and thus the approximate and ex- 
act minima approach zero. We can see a weak correlation between mean and minimum 
values: this happens because both values are strongly influenced by the presence of large 
state chains, but the minimum value depends on lattice details in an even stronger way. 

Large state chains give rise also to a gap between the two lowest-lying peaks. This gap 
develops as -4= approaches an integer from below, giving the A-spectrum a very different 
aspect, as shown in figure [121 As the longest state chain grows, the peak near the origin 
becomes taller and well separated from the second peak. This separation is greater than 
the first peak's width, so that it is effectively isolated from the second peak. 

We can compute an estimate of the gap using the minimum A estimate given in equation 
( I14ip . This estimate will be reliable because the approximations leading to it are valid in 
the two first peaks of the distribution. We can define the gap T as the distance between 
the maximum A of the first peak and the minimum A of the second. Both of them are 
known (see equations (11411) and (I133P ). so we have the following formula for the gap as a 
function of A: 



r(A) = a 



(2nd) 
min 



A 



1st 



2V2 





1 




3 1 




1 




M 




< 




N 


•01 




L^aJ 




~ 2 




L^aJ 





. (146) 



This gap is shown in figure [131 where it is shown with respect to the width of the first peak. 
This width is computed using the standard deviation 01(A) of the o>i(p, A) distribution, 
which is 

, 2, 



a 1 (A) 2 = ((A-(A) Wl )V 

(p- (A) Wl ) 2 wi(p,A)dp 



A/1 (A) 
(0.11A) 2 



(147) 



The first peak of U2 is a u>\ distribution with A replaced by A — (L^J + l) 2 1 an d thus 
its width is given by 

o>t(A) « 0.11 ( A- I — I i ). (118) 
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Average and minimum X versus A (J = 2) 



Average and minimum X versus A (J = 2) 




Figure 11: Comparison between brute- force data (bullets, diamonds) and approximate for- 
mulae (solid, dashed lines) for (A) W2 (bullets, solid lines) and A m i n (diamonds, dashed lines). 
Left panel focuses on a small interval enclosing six integer values of -7=. Valleys of both 
magnitudes at those integer values are caused by long state chains. The wild fluctuation of 
samples is a consequence of the lattice details. The small unevenly-spaced peaks of the solid 
line are there because formula (1145ft lacks a diagonal correction as in figure El Right panel 
shows a much larger interval, but samples have been taken only at half- integer values of -^=. 
The apparently constant profile of (A) W2 is caused by the scale: it is actually decreasing at a 
rate ten times smaller than A m i n . Note how formula (|141h for A m j n (dashed line) works as an 
almost saturated upper bound for A m i n . 




Figure 12: Logarithmic histograms and density of states (solid line) for three near values of 
A characterizing three different examples of the J = 2 Einstein-Maxwell landscape. When 
approaches an integer from below (100 in this case), the first peak shifts to the left and 
becomes isolated, thus creating a gap in the A-spectrum. 
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or some multiple of it. Immediately we can see that this width will approach zero as ^= 
approaches an integer, and therefore the relative gap will become enormous. There is some 
values of A for which the gap becomes negative, that is, the two first peaks of U2 overlap. 
This can happen if -7= > 37, as can be seen in figure [131 This gap will never disappear for 

large values of -7= because the width of the first peak will always vanish at integers, but 
the intervals of positive gaps are smaller when A decreases. That is, the gap is positive 
for ^= < 36, and for greater values the gap changes sign between two consecutive integer 
values of -4=?=, with the zero being closer and closer to [-75=] as A decreases. 




Figure 13: Gap of the J = 2 density of states as a function of A (equation (|146p ). It is 
shown the quotient between the (absolute value of the) gap and the first peak width. Gaps 
can be negative for sufficiently low A, and the negative-gap intervals become greater when 
A becomes smaller. Decreasing curves represent negative values of the gap, while increasing 
ones represent positive values. 
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5 Comparison between the Bousso-Polchinski and multi- 
sphere Einstein- Maxwell landscapes 

Obviously, the multi-sphere Einstein-Maxwell landscape cannot be considered as a model of 
the string theory landscape, because it belongs to a completely different family of theories. 
Nevertheless, the features we have described in the previous sections are not excluded from 
the string theory landscape, and they are qualitatively different in other simplified models, 
such as the Bousso-Polchinski (BP) landscape pp. We will now provide a brief summary 
of the main features of the BP landscape, and then we will stress the differences with the 
multi-sphere EM landscape. 

5.1 The Bousso-Polchinski landscape 

The BP landscape is a simplified model which provides an elegant method for solving the 
cosmological constant problem^. The starting point is M-theory, which is formulated in 
10+1 dimensions, compactified down to 3+1. One of the main ingredients of this theory 
is a seven-form, which is used to introduce the Brown- Teitelboim cosmological constant 
neutralization mechanism [121113], which is a generalization of the Schwinger pair creation 
process responsible for the spontaneous lowering of a strong electric field. 

In the presence of a compactification manifold having three-cycles, the seven-form is 
expanded in a basis of harmonic three-forms, whose coefficients are four-forms. After 
dimensional reduction, the four- dimensional duals of the four-forms are zero-forms, that 
is, scalars, which are quantized by virtue of generalized Dirac quantization conditions. 
The total value of the flux of a four-form in the j th three-cycle is an integer multiple of 
a fundamental charge qj which is proportional to the volume of the three-cycle. These 
charges are moduli of the theory, whose stabilization is given a priori in the BP model. 

A vacuum state of this model is given by specifying the integers representing the value of 
the four-form flux stored in each three-cycle. Transitions between the states are mediated 
by instantons, which can be viewed as M5-brane bubbles with two "legs" enclosing a three- 
dimensional interior of a different vacuum energy density, while having three remaining 
"legs" wrapping the flux in a three-cycle. 

Thus, the vacuum states of the model are arranged in the nodes of a lattice in flux 
space. A given state is specified by J integers n%, ■ • • ,nj, whose effective cosmological 
constant A is given by 



In equation (I149p . J represents the number of three-cycles inside the compactification 
manifold; A is the bare cosmological constant of the theory, which should be negative, so 
that A can reach a small value; and qj are the moduli, that is, the elementary charges of 
the fluxes. 

6 Good reviews of the cosmological constant problem can be found in references [TCTl [TTj 




(149) 
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There is a Minkowski surface in flux space separating AdS and dS states, which is 
obtained by setting A = in (I149p . A node of the lattice can be located very close to 
this surface, and the number of such nodes can be huge by choosing a large enough J. 
The existence of these nodes, randomly close to the A = surface, is essentially the BP 
mechanism solving the cosmological constant problem. 

In a BP landscape with a large amount of fluxes, the vast majority of the nodes are 
located far away from the origin. Some criterion is needed to limit the value of the integers 
rij and render the landscape finite. Usually, this is accomplished by introducing a cut-off 
Acutoff in flux space which characterizes the maximum value of A to be possibly reached. 
The computation of the probability of a given state among all available states based on 
abundance of states gives very small values for a large J, and a large J is needed to reach 
a value of A as low as the observed value 10~ 120 [HJ SS]- So this model has the neces- 
sary states, but a very low probability for them to be occupied, which leads to anthropic 
arguments. 

The anthropic window is an interval of cosmological constant values which allow the 
formation of observed structures (like galaxies, stars and planets) [16]. Even inside this 
anthropic window, the number of states is so huge that the probability of a state having 
A = 10~ 120 is tiny. Dynamical relaxation inside the BP landscape reduces the states 
to a shell wider than the anthropic window [47], and thus do not solve this problem. 
This requires less convincing anthropic arguments to explain the observed value of the 
cosmological constant. 

We can rephrase this problem by saying that the distribution of A values near A = is 
flat [HI HQ] , and thus specially small values do not get rewarded. The observed value of 
A lies in a very thin shell, very small when compared with the anthropic or dynamically 
relaxed shells. Thus, the flatness of the distribution gives rise to such small probabilities. 

In addition, there is another complication with large values of J. When the dimension 
of flux space is large, the vast majority of states in any spherical shell are confined to 
coordinate hyperplanes with a dimension of near J a* with a* < 1 [29]. The bulk of 
the spherical shelfj] is almost devoid of states, and the number of non-vanishing fluxes is 
generically less than J. Nevertheless, stability arguments often force the integers rij to be 
nonzero, even large ones; this would dramatically lower the number of states in the BP 
landscape, resulting in an empty anthropic shell. This a*-problem of the BP landscape is 
not restricted to sets with spherical symmetry; secant states are not spherically distributed, 
and present the same problem. 

5.2 Compasison between BP and ms-EM landscapes 

The previously described features of the BP landscape contrast with their counterparts in 
the multi-sphere Einstein-Maxwell landscape. First of all, this landscape is derived from a 
2 J + 2-dimensional theory, after its dimensional reduction to 1+1 dimensions. Therefore, 

7 That is, the region of the spherical shell surrounding the diagonals of flux space, where states are 
located far away from the coordinate hyperplanes. 
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the resulting cosmologies are not comparable. Nevertheless, we will focus in the distribution 
of states and qualitative features of the landscape. 

The "bare" cosmological constant A is negative in the BP case, allowing cancellation 
in the effective cosmological constant A; if A were positive, no AdS nor low- lying dS states 
would longer exist. In the EM case, A should be positive; otherwise, dS states would 
not exist at all. Thus, both landscapes have twin versions with reversed A which are not 
physically interesting. 

The BP model assumes that its moduli are frozen by some external, unspecified mecha- 
nism. Therefore, the elementary charges are parameters of the model, as well as A. On the 
other hand, the moduli of the EM theory, which are the radii of the internal spheres, are 
fixed (at least at a linear level) by an effective potential built from the magnetic field, the 
curvatures and the vacuum energy density. Thus, this theory needs only one parameter, 
A. It is generally believed that the same stabilization mechanism should work in the BP 
model, but as far as we know it has not been implemented yet. 

The simplicity of the formula for A in the BP model, (I149p . is to be compared with the 
equation determining A in the EM model, f )68|) . (198]) or fll02j) . In this equation, A cannot 
be isolated in general, and there are several branches for each node. Nevertheless, only 
the principal branch has solutions with positive curvatures, and these are the only ones 
with a chance of being stable. Moreover, this equation can have zero, one or two solutions, 
depending on A, giving zero, one or two states per node in flux space. In contrast, (I149p 
always has one solution, and no more, per node. Thus, the correspondence between nodes 
and states is one-to-one in the BP model, but this is not the case in the EM model. 

The finiteness of the BP landscape is a consequence of a cutoff introduced in flux space. 
As commented above, were this cutoff absent, the theory would have an infinite family of 
states with infinitely high-A, which would raise the problem of choosing initial conditions. 
The EM model has a finite amount of stable dS states because of the presence of a branching 
point in the equation determining A. The number of unstable dS states is much greater, 
but these states are excluded from the landscape. Thus, the stability analysis gets rid of 
the majority of dS states, and so we expect this situation to be analogous in the BP model 
completed with a stability analysis. This ingredient can thus significantly change a lot the 
general properties of the BP model, because it would exclude a huge amount of states from 
the landscape. This might be a feature, though, because it might raise the probability 
for the system to be in an anthropic state, which is currently very small because of the 
enormous amount of dS states present. But it is impossible to establish this claim or the 
opposite without a well-defined model to work with. 

AdS states are finite in number in the BP model, because they are located inside a 
sphere in flux space. In the EM model with J > 1 there is an infinite number of them 
and they are always stable. Therefore, the probabilistic arguments based on the number 
of states cannot be applied here, because the probability of dS states would always be 
zero. This argument might be interpreted as indicating that the method of computing 
probabilities using simply amounts of states could be completely wrong in both models. 
As a consequence, the probability measure used in these landscapes should be revised from 
scratch. 
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Both models have a A = surface separating dS from AdS states, which in the BP 
model is 

2|A| = X>X 2 , (150) 

3=1 

that is, a sphere in flux space (parametrized in qjUj coordinates), while in the EM model 
it is 

A = E^> (151) 

which is a sphere after performing a coordinate inversion. These surfaces provide the 
BP mechanism for solving the cosmological constant problem: if the landscape contains 
a state randomly close to this surface then this state can have a realistic value of A. 
Both models have this property. Nevertheless, the surface (115ip is not compact, and it 
allows for long state chains whose cosmological constant can approach very small values. 
This phenomenon is absent in the BP model, and constitutes a basic difference because it 
increases the amount of states in the anthropic shell. As stated above, we don't know the 
correct way of computing probabilities, but state chains provide a new source of low-lying 
states which is absent in the BP model. 

State chains are also responsible for a crucial modification in the distribution of A 
values, which is flat near A = in the BP stated above. In the EM case, the a; (A) 

density vanishes at A = 0, but it has a huge peak of small values, corresponding precisely 
to those lying in the state chains. Thus, this distribution is not flat, which means that the 
randomly-close-state mechanism is less important in the EM model than in the BP case, 
because only states near the diagonal in flux space contribute to it, while all states near 
the sphere contribute in the BP case, accounting for the difference. State chains provide 
a dominant peak of small A values, which is another different mechanism for solving the 
cosmological constant problem. This mechanism has its own drawback, because this peak 
can be very near A = for very specific values of A (namely, when -4= is very close to an 
integer from below), but this is not the case for generic A values. Thus, both mechanisms 
are different, and both can be fine-tuned to have very small values of the cosmological 
constant. Quantifying the neccessary fine-tuning is not easy without a valid method for 
computing probabilities. 

Finally, the a*-problem is absent in the EM model, because the A = surface never 
approaches the coordinate hyperplanes where one or more rij = 0. Thus, if the stability 
results found in the EM model translate to the BP model completed with a stability 
analysis, then we are forced to conclude that the vast majority of dS states, which are 
near the hyperplanes, would be unstable, and thus there would be excluded from the BP 
landscape. This would change all reasoning based on number of states, if it were to be of 
any use, and thus making further anthropic assumptions would be totally unneccessary. 

Table [1] summarizes all the issues we have addressed while comparing the BP and EM 
landscapes. 
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BoussoPolchinski 



Multi-sphere Einstein-Maxwell 



Parameters 



A > Uj}j=i -J 



Sign of A 



A < 



A > 



Stability 



Assumed 
(qj given a priori) 



Linear stability 
understood (k > 0) 



A formula 



2 2 



A = L n (X) (principal branch) 



j=i 



Finiteness 



A < A 




A < A b (dS) 
Infinite (AdS) 



A = surface 



J 1 



Cosmological Randomly 
constant problem close states 



Randomly 
close states 

and 
Long state 
chains 



l«2 






111 



cj(A) distribution 




Dominant 
and many 
spikes 




a* problem 



Vast majority of states have 
(1 — a*) J vanishing fluxes 
(and they are probably 
unstable) 



There are no stable states with 
vanishing fluxes 



Table 1: Summary of the differences between the Bousso-Polchinski and multi-sphere 
Einstein-Maxwell landscapes. The first two graphics emphasize the different distribution of 
states in flux space, which is the origin of the state chains. These state chains are responsible 
of the dominant spike in the A-density, as shown in the two last graphics. 



51 



6 Possible implications for the string theory land- 
scape 

The comparison carried out in the previous section leads to some features that a complete 
treatment of the BP landscape interpreted as a toy model of the true string theory land- 
scape would bear when compared to what is currently believed. This section is devoted to 
envision what the BP landscape would look like if some of the main features of the EM 
landscape were found to hold. 

Two key points should be stressed: 

• Theories with four-form fluxes have duals which are gauge-gravity theories. The 
main difference between them is that gauge-gravity theories have a built-in cut- 
off mechanism which limits the maximum values of the integers characterizing the 
landscape. Thus, those landscapes have natural finiteness conditions, such as the 
branching point in the multi-sphere Einstein-Maxwell landscape. This would avoid 
the neccessity of a cut-off scale put by hand in the model. 

• Stability conditions are important not only beacuse they complete the model, but also 
because they exclude a huge amount of states from the model. If the same exclusion 
were to take place in the BP landscape, the vast majority of its nodes would not 
be true states of the landscape, and the predictions of existence of anthropic states 
would change dramatically. Thus, stability conditions have a two-fold purpose: on 
the one hand they fix the values of the elementary charges to be used in the model, 
and on the other hand they limit which nodes have physically relevant states. A 
priori frozen moduli fulfill the first purpose, but do not help in deciding which nodes 
have states, which causes a huge proliferation of states, which may be spurious. The 
EM model shows that the vast majority of dS states are unstable, and so we can 
expect the same to be true in a completed BP landscape. 

Other features of the EM landscape may not have a direct translation to a completed BP 
model, such as state chains. They are a consequence of the asymptotes found in the null-A 
surface, which is non-compact. This is an indication that the details of the null-A surface 
provide different sources of low-lying states which change the density of states a; (A). This 
distribution is needed when one has to compute probabilities in a given landscape model; 
but a completely clear, unambiguous, quantum prescription for computing probabilities 
with a general model is still lacking, and therefore the implications of the details of cu(A) 
in the computation of probabilities cannot go beyond the naive arguments based on state 
abundances. At this simple level, details of the null-A surface translate in peaks in the oj(X) 
density, thus producing very different probabilities for the states in the anthropic shell. The 
BP model has a spherical null-A surface, and thus ou(X) has no peaks; if the null-A surface 
of more realistic flux compactifications of M-theory had other nontrivial shapes, this would 
be reflected in the a; (A) distribution and in the final computation of the probabilities. So 
this is the last point that the EM landscape brings in: the details of the null-A surface are 
very important for probability computations. 
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7 Conclusions 



We have addressed a simple sector of the Einstein-Maxwell theory as an exactly solvable 
model of a landscape. The theory, formulated in 2J + 2 spacetime dimensions, has a 
single parameter in the Lagrangian, namely, the "bare" cosmological constant A > 0. 
The compactification has the form (A)dS 2 x (S 2 ) J , which is referred to as multi-sphere 
Einstein-Maxwell compactification. Equations of motion for the corresponding metric 
ansatz are algebraic equations for the values of the curvatures of the inner spheres and the 
effective cosmological constant A of the cosmological part. In the presence of a magnetic 
monopole, the magnetic flux in each sphere, which is quantized by a Dirac condition, 
stabilizes the configuration which spontaneously would decompactify. The cosmological 
constant A helps to evade the Maldacena-Nunez no-go theorem in this case [^9l 123] . The 
different combinations of the flux quanta stored in the spheres give rise to a complicated 
landscape, in which each configuration of integers (called a node) can host a true stable 
state of the model, two stable states, an unstable state and a stable one, or no state, giving 
rise to two branches of (AdS and dS) states. It is found that for J > 1 an infinite family of 
stable AdS states exist, but stable dS states exist only near the branching surface, which is 
the locus at which both branches meet, that is, the null- A surface. The structure of the null- 
A surface gives rise to the state chains, which provide a different source of low-lying states 
besides the randomly close states which help to solve the cosmological constant problem in 
the Bousso-Polchinski landscape. State chains also help in counting states approximately, 
and they translate in peaks in the density of states cu(A). 

All the previous features of the model are qualitatively different from its counterparts 
in the Bousso-Polchinski landscape. We think that in a completed BP model, all these 
differences would render a very different picture with respect to the number of states, 
probabilities and anthropic reasoning. Thus, despite not being a realistic landscape model, 
the multi-sphere Einstein-Maxwell model has very appealing features that might propagate 
in more realistic models of the true string theory landscape. 

A fundamental piece missing in this multi-sphere Einstein-Maxwell landscape is the 
cosmological measure. As a very simple 1+1 cosmological model, it is plausible that this 
problem would be easier to solve than its 3+1 counterpart in more complicated compacti- 
fications. This problem will be addressed in future papers. 
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A Effect of derivative couplings in the multi-radion 
evolution equations 

In this appendix we give a heuristic argument leading to the conclusion that the linear 
stability analysis of the multi-radion field evolution equations, equation (185]) . which is 
achieved by neglecting the derivative couplings, can be promoted to a non-linear stability 
analysis in which the linear stability is preserved as long as perturbation amplitudes are 
sufficiently small. 

To begin with, we consider again equation ( 1851) : 



k k 



X- e -^^u'AQ. (152) 



The derivative couplings appear in a quadratic form. We will use the symbol £ to denote the 
J-component column vector of the perturbations £j, and then we will write the derivative 
couplings in matrix form as 



- e- 2 V" 



(153) 



where the constant J x J matrix Mj has the number —1 along the diagonal except for 1 
along the j th row and j th column and zeroes elsewhere: 
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(154) 



The matrix Mj has a spectrum with two simple eigenvalues ±v J and a (J — 2) -degenerate 
— 1 eigenvalue. It is therefore an "almost negative-definite" quadratic form, which repre- 
sents a non-linear force. 

We are interested in the effect of the derivative couplings; therefore, we approximate 
the system of equations by Taylor-expanding to first order the right-hand side of equation 
( I153p . thus considering the much simpler system (written in matrix form) 



- H jk £ k 



(155) 



where the frequency matrix H is given by equation 0881) in section |X 
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As said above, a complete analysis of the system f II 5 5 j) is not possible, but we can form 
a single equation out of it by projecting the system of equations along a constant direction 
w = (wj): 



(w,£>a0 + (£«,M w £g) 



The matrix M w is 



(w,M) = £> 3 -M,-. 



(156) 



(157) 



Any solution of the system (I155p is a solution of equation (11561) . but the reciprocal is false. 
Nevertheless, we can look for normal modes of the form £(t, x) = x{t, x ) v f° r constant v 
and a single scalar function x(t, x)\ 



- e- 2 ^ 



(w, v)Xa/3 + (V, M w v)XaXl3 



(w,Hv)x- 



Taking x — x{t) an d w = v for simplicity, we obtain 
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(158) 



(159) 



If the non-linear term were absent, we would have a simple oscillator equation with a 
frequency given by the Rayleigh quotient of the matrix H. The solution of this equation 
would be a solution of the linear system if v were chosen as an eigenmode of H. In this 
sense, the projected equation (II 59ft is an average equation, and its solution (a weak solution 
henceforth) can indicate the behavior of the true solutions we are inspecting. Of course, 
this is a heuristic argument, but we can argue that true solutions provide weak solutions; 
thus, an unstable true solution should be reflected by an unstable weak solution. This 
argument has the obstacle of the existence of the normal modes we are using as ansatz; 
thus, as long as the normal modes constitute a reasonable description of the system 



the projected equation will reflect accurately the character of its non-linear counterpart. 

Nevertheless, we can use the projected equation (I159P to see if the presence of the 
non-linear term can render unstable a linearly stable solution. 

We will rewrite equation fl!59[) as 

{V ; M » V) , and h = e*> ^ gff ■ 
(v,v) {v,v} 



X + rnx 2 = —hx , with m 



(160) 



The parameter m depends on the projection direction v but not on time. In contrast, h 
depends also on time by the presence of the e 2 ^ factor. We will now discuss the expected 
domain of both parameters in the following. 

• The parameter m is the Rayleigh quotient of the matrix M v on the projection direc- 
tion v. The eigenvalues of M v are ±y/J v 2 , both of them nondegenerate, and 

— . Vj with J — 2 degeneracy. Thus, depending on the projection direction v, m 
can have both signs. If v is taken to have unit norm, then m will be some value in 
the interval [— vJ]. 
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• The 1+1 cosmological solution <p(t) has a characteristic evolution time t<f> = ^p^- O n 

the other hand, the longest characteristic evolution time of the oscillator term is = 
-^=, where k is the minimum eigenvalue of the frequency matrix H. Of course, we are 
considering a linearly stable equation, so that k > 0. We can consider as slowly- 
varying if its characteristic time is much greater than the oscillator characteristic 
time: 

fy»*H =^ A<k (161) 

Thus, for times t ~ we may consider as constant, and thus will be a positive 
number. 

Equation (11601) is integrable. We can show the form of its trajectories by writing \ = 7: 

* = 7 1 d ^ r X / 1RO x 

2 > — = —my — h— . (162) 

7 = ~ m 7 - ^Xj dx 7 

The last equation has the exact solution 



^ =^-s?J e "^ + mfe-*;- (l63) 

In the previous equation, 70 = 7(0). Specializing m = (which removes the nonlinear 
term) we obtain the oscillator trajectory 7 2 + hx 2 = 7q- The points where 7 = are called 
turning points of the trajectory, and they mark its domain because of the square in (11631) . 
The trajectory has one or two turning points given by the equation 

v2_ J_V*»x = Af v - J-^ 



^° 2m 2 ) m\ 2m J ^ ^ 

If m is a fixed positive value, then the previous equation has a single solution if 7q — 
2^2 > but it has two solutions if y 2 — ^2 < 0. A single turning point describes an 
open trajectory, while two turning points describe a closed one, see figure TAL. Thus, the 
oscillator trajectory remains closed when we turn on the non-linearity if 

* < 2^ (165) 

Therefore, a linearly stable trajectory remains non-linearly stable if the amplitude 70 does 
not exceed the critical value = • Beyond this value, the trajectory is open and 
therefore the linearly stable solution becomes non-linearly unstable. 

The lowest value of the critical amplitude is reached when m is largest; for a unit-norm 
projection direction, the largest value of m is yJ, as discussed above. Thus, the lowest 
value of the critical amplitude is 

7c 2 min = ■ (166) 
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Turning points of the x, V trajectory Trajectories of the pseudo-normal mode in the %, Y plane 




Figure 14: Turning points (left) and actual trajectories (right) of the normal modes of the 
projected multi-radion evolution equation. 



As long as amplitudes are smaller than this value, the non-linear derivative couplings 
cannot spoil linear stability. Nevertheless, when k is small, perturbations have a chance of 
trigger a non-linear instability and destabilize a linearly stable state. 

Summarizing, the non-linear terms in the multi-radion evolution equations respect the 
linear stability criterion except in the regime of large amplitudes, which is most easily 
accessible when the minimum eigenvalue of the frequency matrix becomes small, that is, 
in the onset of instability. 
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